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ABSTRACT 

We present an analysis of XMM-Newton observations of the superbubble 30 Dor 
C and compare the results with the predictions from the standard wind-blown bubble 
model. We find that the observed X-ray spectra cannot be fitted satisfactorily with 
the model alone and that there is evidence for nonthermal X-ray emission, which is 
particularly important at ^ 4 ke V. The combination of the bubble model and a power- 
law gives a reasonable fit to the observed spectra. The thermal pressure and central 
temperature of the bubble are 3.3 x 10~ n dyne cm -2 and 7.4 x 10 6 K, respectively, 
and we infer that, for a bubble age of t ~ 4 x 10 6 years, the ambient density is no — 
38 cm -3 , the mechanical luminosity is £ m cch ~ 10 37 erg s , and the coefficient of 
thermal conductivity is ~ 0.05 of the Spitzer value. The total unabsorbed 0.1-10 keV 
luminosities of the eastern and western parts of the bubble are ~ 3 x 10 36 erg s _1 and 
~ 5 x 10 36 erg s _1 , respectively. The unabsorbed 0.1-10 keV luminosity of the bubble 
model is ~ 4 x 10 36 erg s _1 and so the power-law component contributes between 1/3 
and 1/2 to the total unabsorbed luminosity in this energy band. The nature of the 
hard nonthermal emission is not clear, although recent supernovae in the bubble may 
be responsible. We expect that about one or two core-collapse supernovae could have 
occured and are required to explain the enrichment of the hot gas, as evidenced by 
the overabundance of a-elements by a factor of > 3, compared to the mean value of 
~ 0.5 solar for the interstellar medium in the Large Magellanic Cloud. As in previous 
studies of various superbubbles, the amount of energy currently present in 30 Dor C is 
significantly less than the expected energy input from the enclosed massive stars over 
their lifetime. We speculate that a substantial fraction of the input energy may be 
radiated in far-infrared by dust grains, which are mixed with the hot gas because of the 
thermal conduction and/or dynamic mixing. 

Subject headings: acceleration of particles — stars: winds, outflows — ISM: bubbles - 
HII regions — X-rays: individual (30 Dor C = DEM L 263) — X-rays: ISM 

1. INTRODUCTION 

The fast winds from massive stars and their subsequent supernovae (SNe) ejecta dominate the 
mass and energy input to the interstellar medium (ISM; e.g., Abbott 1982) and so may regulate star 
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formation and galaxy evolution. Most massive stars form in OB associations, and the concentrated 
energy input from massive stars in OB associations sweeps up the ambient ISM to form expanding 
shells called superbubbles (e.g., McCray & Snow 1979). The physical structure of a superbubble is 
very similar to that of a bubble blown by the stellar wind of an isolated massive star (Mac Low &; 
McCray 1988). However, theories of such bubbles are still very uncertain and have not had much 
serious appraisal with observations of high quality. 

Studies of Galactic superbubbles are difficult because of their large angular sizes and the 
confusion of objects in the line of sight. The Large Magellanic Cloud (LMC), in contrast, provides 
an excellent laboratory to study superbubbles since it is nearby (at a distance of 50 kpc, as adopted 
here, 1" = 0.24 pc; Eastman & Kirshner 1989) and almost face-on (inclined by ~ 35° to the line of 
sight; van der Marel & Cioni 2001). It also has little foreground (E-Q-y = 0.08 mag; Xu et al. 1992) 
or internal extinction (£b-v m the range from ~ 0.06 to ~ 0.3 mag; Bessel 1991). The superbubble 
30 Dor C (= DEM L 263; Davies, Elliott, & Meaburn 1976) is southwest of the main 30 Doradus 
complex in the LMC, and is coincident with the large OB association LH 90 (= NGC 2044; Lucke 
& Hodge 1970). This OB association consists of a number of stellar clusters and is particularly rich 
in Wolf-Rayet (W-R) stars (Lortet & Testor 1984). The most massive stars appear to have ages in 
the range from ~ 3 to ~ 7 Myr (Testor, Schild, & Lortet 1993). In X-ray imaging observations, 30 
Dor C appears almost as a complete shell with a diameter of ~ 7' (Wang & Helfand 1991; Dennerl 
et al. 2001; Dunne, Points, & Chu 2001). The X-ray shell is confined within the Ha-emitting 
filaments (see Fig. 11 of Dunne et al. 2001), as expected from the evaporation of cool gas into the 
hot interior of a superbubble. There is also an indication for a substantial nonthermal contribution 
in addition to an assumed one-temperature thermal component, particularly in the western part of 
shell (Bamba et al. 2003, 2004). The radio emission from 30 Dor C may be due to a combination of 
synchrotron-emitting relativistic electrons produced by shocks within the superbubble and thermal 
emission from the Ha-emitting shell (Mathewson et al. 1985). 

In this paper, we report on a spatially resolved spectroscopic analysis with XMM-Newton of 
the diffuse X-ray emission emanating from 30 Dor C. We specifically use the wind-blown bubble 
model of Weaver et al. (1977; see also Mac Low Sz McCray 1988) to infer the temperature and 
density structure of the X-ray emitting gas, which in turn are used to estimate both the density of 
the ambient medium and the rate at which energy is injected into the ISM by the stellar winds. 

2. OBSERVATIONS AND DATA REDUCTION 

The region surrounding 30 Dor C is known to be rich in X-ray sources and was chosen as 
the first light image for XMM-Newton (Jansen et al. 2001). However, data acquired this early in 
the mission are in a format which cannot be processed by the Science Survey Center using the 
standard software and so are not included in the XMM-Newton archive at the time of writing. The 
data reported here were instead acquired during ~ 21.6 and ~ 37.5 ks observations of SN 1987A 
and PSR J0537-6909 with XMM-Newton on 2000 November 25-26 (observation ID 0104660301; 
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PI: Watson) and 2001 November 19-20 (observation ID 0113020201; PI: Aschenbach), respectively. 
These observations had a low background count rate and were pointed ~ 6' (0104660301) and 
~ 9' (0113020201) offset from the center of the 30 Dor C. While there are other observations that 
covered 30 Dor C in the XMM-Newton archive, these observations are either severely contaminated 
by background flares (observational ID 0083250101) or had a very short good exposure time (e.g., 
~ 3.3 ks for observation 0104660101; PI: Watson). The EPIC MOS cameras (Turner et al. 2001) 
were operated in the full- frame mode with the medium filter inserted. Over the time interval 
between the two observations, the energy resolution of MOS CCDs had worsened, from 140 eV to 
150 eV (FWHM) at 5.9 keV, due to radiation damage. These values are actually worse than the 
initial post-launch value of ~ 130 eV. The on-axis telescope point spread function is 4."3 (FWHM) 
at 1.5 keV and the field-of-view is roughly circular with a diameter of ~ 30'. Unfortunately, the 
EPIC pn exposure of the first observation was taken with the calibration lamps switched on, so 
rendering the data useless for scientific analysis, and the pn was operated in the timing mode during 
the second observation, so providing no imaging data. 

The data reduction and analysis were done using the XMM-Newton Science Analysis Software 
(SAS) version 5.4.1 (released on 2003 January 16). The MOS data were reprocessed into calibrated 
event list files using the SAS program emchain. The procedure of Read & Ponman (2003) was used 
to filter each event list for periods of high background. We created lightcurves in the 10-15 keV 
energy range, where no source counts are expected due to the very low effective area of the X-ray 
telescope, and using only single pixel events and flag values as defined by #XMMEA_EM. Periods 
when the count rate in these lightcurves exceeded 0.35 counts s _1 were excluded from the event 
files. The event lists were then filtered further, keeping only single, double, triple, and quadruple 
pixel events in the 0.2-15 keV energy range. For the first observation, the total good exposure time 
was ~ 21 ks, and for the second observation, the MOS1 and MOS2 CCDs had total good exposure 
times of ~ 17 ks and ~ 13.5 ks, respectively. We also checked for soft proton flares in the 2-8 keV 
band, but found none. 

The interstellar extinction (£b-v) towards the western part of 30 Dor C is known to be higher 
by at least 0.4 mag than that towards the eastern part (Testor et al. 1993; see also Dunne et al. 
2001 and Fig. 5 of Dennerl et al. 2001). We therefore extracted spectra from two sectors of an 
ellipse centered on 30 Dor C with major-axis diameter of ~ 6.'9 (~ 100 pc) and major-to-minor 
axis ratio of 0.9. The regions are marked in Figure 1, which is an image of the central region of the 
30 Dor C field (see § 3.1 for a description of X-ray sources in this image). The shape and size of 
the extraction regions were chosen to include as much of the diffuse emission from the superbubblc 
as possible, while minimizing the contribution from the background. 

The background was estimated from source free regions near the superbubble (Fig. 1). There 
are spatial variations in the detector background rate, which can be estimated from a comparison 
of blank-field spectra in the source and source free regions of the 30 Dor C observations. For this 
comparison, we have used blank-field observations with a total exposure of ~ 488 ks and ~ 593 ks 
for the MOS1 and MOS2 detectors, respectively (Read & Ponman 2003). These data, taken with 
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the same instrument mode and filter combination as the 30 Dor C observations, had point sources 
of emission excised, and were reprojected on the sky to match the positions of the 30 Dor C 
observations. 

We find that for the first observation, the MOS1 blank- field spectra in the source free region 
and in the eastern part of 30 Dor C are statistically consistent (within 3a deviation) with each 
other in the 0.3-10 keV band, with x 2 = 573.5 (497 dof). However, for the same observation, a 
comparison of the MOS2 blank-field spectra in the source free region and in the eastern part of 30 
Dor C revealed significant differences in the background count rate near the energy of the Si Ka 
(1.7 keV) fluorescence line. We therefore compared the MOS2 blank-field spectra in the 0.3-1.6 
and 1.85-10 keV band, and find that the spectra are consistent with each other, with x 2 = 510.4 
(466 dof). There are also, for the first observation, and near the energies of the Al Ka (1.5 keV) 
and the Si Ka fluorescence lines, significant differences in the background count rate between the 
blank-field spectra in the source free region and the blank-field spectra in the western part of 30 
Dor C. In the 0.3-1.3 and 1.85-10 keV band, however, the blank-field spectra in the source free 
region are statistically consistent with the blank-field spectra in the western part of 30 Dor C, with 
X 2 = 465.6 (449 dof) and X 2 = 458.6 (430 dof) for the MOS1 and MOS2 detectors, respectively. 

Similarly, for the second observation, and near the energy of the Si Ka fluorescence line, we 
find that there are significant differences in the background count rate near between the blank- 
field spectra in the source free region and the blank-field spectra in the eastern part of 30 Dor 
C. In the 0.3-1.6 and 1.85-10 keV band, however, the blank-field spectra are consistent with each 
other, with X 2 = 501.4 (464 dof) and X 2 = 471.0 (476 dof) for the MOS1 and MOS2 detectors, 
respectively. There are also, for the same observation, and near the energies of the Al Ka and Si 
Ka fluorescence lines, significant differences in the background count rate between the blank-field 
spectra in the source free region and the blank-field spectra in the western part of 30 Dor C. We 
therefore compared the blank-field spectra in the 0.3-1.3 and 1.85-10 keV band, and find that the 
MOS1 and MOS2 blank-field spectra in the source free region should be scaled by 0.92 and 0.84, 
respectively, in order to match the blank-field spectra in the western part of 30 Dor C. The rescaled 
spectra are statistically consistent with the blank- field spectra in the western part of 30 Dor C, with 
X 2 = 546.3 (478 dof) and X 2 = 451.1 (468 dof) for the MOS1 and MOS2 detectors, respectively. 

We therefore excluded data in the 1.6-1.85 keV band from the MOS2 spectra of the diffuse 
emission in the eastern part of 30 Dor C and data in the 1.3-1.85 keV band from the spectra of 
the diffuse emission in the western part of 30 Dor C. For the second observation, we also excluded 
data in the 1.6-1.85 from the MOS1 spectrum of the diffuse emission in the eastern part of 30 Dor 
C. 

We created spectral redistribution matrices and ancilliary response files for the source spectra 
using the SAS programs rmfgen and arfgen, respectively. The ancilliary response file appropriate for 
each source spectrum was created weighting the position-dependent quantum efficiency and effective 
area by the number of 0.3-10 keV blank-field background subtracted counts in the corresponding 
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region. Prior to performing the spectral analysis with XSPEC version 11.2.0au (Arnaud 1996), we 
rebinned the data so that the background subtracted signal-to-noise ratio in each bin is greater than 
4. When performing the spectral analysis, we excluded data with energies greater than 10 keV, 
since the background dominates the total emission from 30 Dor C at energies greater than this. We 
also excluded data with energies less than 0.3 keV, because of the uncertainties in the calibration 
of the instruments. 

3. ANALYSIS AND RESULTS 

3.1. Morphological Properties 

Figure 1 shows the XMM-Newton image of the ~ 26' x 26' (370 pc x 370 pc) region centered on 
30 Dor C in the 0.3-8 keV band. In this image, 30 Dor C is a roughly elliptical feature with major 
axis ~ 7' (100 pc). This feature appears to be the observational manifestation of a bubble blown in 
the ISM by the winds from massive stars in the bubble's interior (Mac Low Sz McCray 1988). Other 
sources of X-ray emission in Figure 1 include SN 1987A, which is located ~ 6' southwest of the 
bubble center and close to a region of diffuse X-ray emission associated with the Honeycomb nebula 
(Meaburn et al. 1993; Chu et al. 1995b). There is also a compact X-ray source, J0536. 9-6913, to 
the east of 30 Dor C, which has a radio counterpart and is probably a background AGN (Haberl 
et al. 2001). The bright source at the left of the image is the Crab-like supernova remnant N157B 
(Wang & Gotthelf 1998). 

In Figure 2, we compare an Ha emission-line image from University of Michigan/CTIO Mag- 
ellanic Cloud Emission-line Survey (MCELS; Smith et al. 1998) to an adaptively smoothed image 
of the inner 15' x 15' region of the 30 Dor C field in the 0.3-8 keV band. The edge-brightened 
X-ray emission is presumably due to an increase in gas density around the edge of the bubble, 
and is spatially enclosed by the Ha-emitting filaments. It is conceivable that the increase in gas 
density is due to the evaporation of cool (~ 10 4 K) gas into the hot interior from the surrounding 
Ha-emitting shell (Weaver et al. 1977). 

We also show adaptively smoothed images of the central 12.'5 x 12.'5 region of the 30 Dor C field 
in the soft (0.3-1 keV), medium (1-2 keV), and hard (2-8 keV) bands (Fig. 3). In these images, 
there is evidence for spectral variations in the X-ray emission across the face of the bubble. While 
the western part of the bubble is easily seen in the medium and hard band images, it is almost 
invisible in the soft band image, clearly due to the X-ray absorption by foreground cool gas (see also 
Dennerl et al. 2001 and Dunne et al. 2001). There is a region of localized CO emission projected 
onto the western part of 30 Dor C (Johansson et al. 1998). Assuming a 12 CO(1-0) intensity of 9.1 
K km s _1 for the cloud (see Table 3 of Johansson et al. 1998) and a CO to H2 conversion factor of 
1.3 x 10 21 molcm~ 2 (Kkms -1 ) -1 for the LMC (Israel 1997), we estimate an H2 column density of 
A^h 2 — 1-2 x 10 22 molcm~ 2 towards the western part of 30 Dor C. This absorption column density 
is sufficient to absorb the soft X-ray emission from the western part of bubble. Intriguingly, in the 
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hard band image, the western part of the bubble appears much brighter than the eastern part, and 
is presumably more luminous. 

Finally, we note that J0536. 9-6913 is almost invisible in the soft band image, due to either 
absorption by gas in the LMC or absorption intrinsic to the source. This source also went undetected 
in the PSPC observations reported by Dunne et al. (2001), which would require either an increase 
in the absorption column density towards the source or a decrease in the intensity of the source. 



In earlier studies of 30 Dor C with Einstein and ROSAT, the X-ray luminosity expected from 
the stellar wind-blown bubble model of Weaver et al. (1977) was often compared with that derived 
from spectral fits to the data using a single temperature thermal plasma model (e.g., Dunne et al. 
2001). The main disadvantage of this approach is that the data are not directly compared to the 
model and so any test of the model will depend on the measurements used to derive the expected 
X-ray luminosity of the bubble (e.g., the number density and expansion velocity of the H n region). 
Also, the wide range of gas temperatures inside the bubble gives rise to soft X-ray emission in excess 
of that expected from a single temperature plasma model and emission lines from a wider range of 
ionization states than those expected from a single temperature thermal plasma model, which will 
naturally result in an underestimate of the heavy element abundances. The shortfall in the soft 
X-ray emission below ~ 0.7 keV from the single temperature thermal plasma model can also lead to 
an underestimate of the absorption column density, and hence the soft X-ray luminosity. In a more 
sophisticated analysis of 30 Dor C with the Einstein data, Wang & Helfand (1991) compared the 
IPC spectrum with that expected from integrating the model over the bubble interior and folding 
the result through the IPC response. However, with the limited quality of the IPC data, they 
were unable to constrain tightly the basic spectral parameters such as temperature and absorption 
column density, let alone to test other interesting parameters such as the efficiency of thermal 
conduction and metal abundance of the X-ray emitting gas. In a more recent analysis of the 
Chandra and XMM-Newton data, Bamba et al. (2003, 2004) discovered nonthermal X-ray emission 
from 30 Dor C, although the discovery was made using a single temperature, non-equilibrium, 
ionized plasma model to describe the thermal emission from the bubble. 

Here we have compared the X-ray spectra of 30 Dor C directly with the stellar wind-blown 
bubble model of Weaver et al. (1977; see also Mac Low & McCray 1988) by implementing it into the 
spectral analysis package XSPEC. The X-ray emission in their model arises primarily from a region 
of thermally evaporated material (see Fig. 3 of Weaver et al. 1977). The temperature and density of 
the hot gas as a function of radius r within this region can be approximated as T(r) = T c (1— r/i?) 2 / 5 
and n(r) = n c (1 — r/i?) -2 / 5 , respectively, where T c and n c are the corresponding central values, 
and R is the outer radius of the bubble. The X-ray flux, Fx, expected from the bubble is 



3.2. Implementation of the Wind-Blown Bubble Spectral Model 




(1) 
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where Ax(T) is the mean emissivity of the hot gas in the X-rays and D is the distance to the source. 

We have calculated the X-ray spectrum with a multi-temperature, variable heavy element 
abundance MEKAL (Mewe, Kaastra, & Liedahl 1995) model, as implemented in XSPEC, using 

dEM(T) = ^fj^ = K [2(T/T C ) 2 - {T/T,)- 1 ' 2 - (T/T c ) 9 / 2 ] dT/T c (2) 

for the differential emission measure dEM(T)/dT, which is a function of temperature. K [= 
107r.R 3 n 2 /(47rZ) 2 ) = 15Vnl/(8irD 2 )] is the normalization of the spectrum in XSPEC. In this form, 
the model spectrum can be applied approximately to a slightly aspherical bubble of an equivalent 
volume V. 



3.3. Spectral Analysis 

The MOS1 and MOS2 spectra of the diffuse emission in the eastern and western parts of 30 
Dor C (a total of 8 spectra) were simultaneously fitted to the above wind-blown bubble model 
and a foreground absorption of a column density of Ah- The absorption cross-sections and atomic 
abundances were taken from Morrison & McCammon (1983) and Anders & Grevesse (1989), re- 
spectively. While the foreground absorption was allowed to be different between the eastern and 
western parts of 30 Dor C, we constrained the bubble model to have the same temperature, metal 
abundance, and intensity in the eastern and western parts of 30 Dor C. This model, which has 5 
free parameters, gives a poor, simultaneous fit to the MOS1 and MOS2 data with x 2 = 842.8 (349 
dof). We also did not obtain a significant improvement in the fit when we allowed the normalization 
of the bubble model to be different between the MOS1 and MOS2 datasets. 

An inspection of the spectra indicates that much of the deviation from the best-fit bubble 
model is due to an excess of hard X-rays in the spectra of the diffuse emission, particularly in the 
western part of 30 Dor C. Such an hard X-ray excess can naturally be explained as nonthermal 
X-ray emission, which can typically be described by a power-law. We therefore added a power-law 
component of photon index T to our model fits. The normalization of the power-law continuum was 
allowed to vary between the eastern and western parts of 30 Dor C, to account for the spatial vari- 
ation in the hard X-ray emission (Fig 3). In our initial modeling of the data, the metal abundance 
was poorly constrained, with a lower limit of ~ 0.1 solar and an upper limit of several times solar. 
We therefore fixed the metal abundance at the known mean value of 0.5 solar for the ISM in the 
LMC (Russell & Dopita 1992). In comparison to the bubble model, this model (model B+P) gives a 
much improved, although still statistically unacceptable, simultaneous fit to the MOS1 and MOS2 
data with \ 2 = 459.8 (347 dof; Table 1, column 2). The absorbed 0.5-10 keV flux in the western 
and eastern parts of 30 Dor C are ~ 1.2 x 10~ 12 and ~ 1.3 x 10" 12 erg s 1 cm 2 , respectively. 
These fluxes are slightly larger than the fluxes reported in Bamba et al. (2003) for their regions 1-4, 
due to the differences in the sizes of the regions used for extracting spectra between the present 
work and Bamba et al. (2003). We note that a marginal (at greater than ~ 98% confidence — a 2a 
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result) improvement in the fit is obtained when we allowed the normalization of the bubble model 
to be different between the eastern and western parts of 30 Dor C. 

The spectra of the diffuse X-ray emission in wind-blown bubbles have often been fitted to 
single temperature plasma models in collisional ionization equilibrium (e.g., Dunne et al. 2001), 
although with little physical justification. In order to compare with these earlier observations and 
analyses, we have replaced the bubble model with a single temperature MEKAL model. We have 
also fixed the metal abundance of the MEKAL at 0.5 solar, since it is poorly constrained by the 
data. This model (model M+P) gives a similar quality fit to the data as the B+P model, with 
X 2 = 463.7 (347 dof; Table 1, column 3). The best-fit equivalent hydrogen column densities and 
normalizations of the power-law continua are comparable to the values obtained from the fitting 
the data to the B+P model, and the best-fit temperature of the MEKAL is consistent with the 
value reported in Bamba et al. (2003) for the thermal emission in their region 1. 

We have also considered the possibility that the hard X-ray excess in 30 Dor C is attributable 
to a second thermal component that might be caused by a recent supernova explosion, for example. 
We replaced the power-law component in the above M+P model with a second MEKAL model, 
with the metal abundances of the two MEKALs fixed at the same value of 0.5 solar and the 
normalization of the high-temperature MEKAL allowed to vary between the eastern and western 
parts of 30 Dor C. We find that this model provides a much worse fit to the data than do either the 
B+P or M+P models, with x 2 = 550.7 (347 dof). This is due to the data at high energies being 
much too flat compared to the model. Although a significant (at greater than 99% confidence) 
improvement in the fit is obtained when we allowed the metal abundances of the two MEKALs to 
vary together, this model, with one extra free parameter, still provides a worse fit to the data than 
do either the B+P or M+P models, with \ 2 = 509.2 (346 dof). Moreover, the metal abundance of 
~ 0.04 solar is uncomfortably small for the ISM in the LMC. More sophisticated models, such as 
the non-equilibrium, ionized (NEI) plasma codes as implemented in XSPEC (see e.g., Borkowski, 
Lyerly, & Reynolds 2001 for a detailed critique of these models), also gave poorer fits to the data 
than do either the B+P or M+P models. We therefore conclude that the hard X-ray excess in 30 
Dor C is best described by a nonthermal emission component. 

Because 30 Dor C is a massive star forming region with possible recent core-collapse supernovae 
occurring inside the bubble, we have further examined the data for an over-abundance of the a- 
process elements compared to Fe-like elements. We did this by dividing the elements into two 
groups: the a-process elements O, Ne, Mg, Si, S, Ar, and Ca; and the remaining Fe-like elements 
C, N, Na, Al, Fe, and Ni. Elements in each group were constrained to have the same abundance 
relative to solar. The abundances of He and the Fe-like elements were fixed at the solar and 0.5 
solar values, respectively, which are appropriate for the ISM in the LMC (Russell & Dopita 1992). 
This model (model vB+P) further improves the fit to the MOS1 and MOS2 data with x 2 = 430.6 
(346 dof; Table 1, column 3). For clarity, we show the co-added MOS spectra from the eastern 
and western parts of the bubble, together with the residuals to this model fit (Fig. 4). There are 
some large residuals around the energies of the Al Ka fluorescence line that arises from the detector 
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housing and the Au M-edge that arises from X-ray absorption in the telescope mirror surface. With 
the possible exception of Si, the a-process elemental abundance of (1.6^05) solar (errors here and 
elsewhere in this paper are 90% confidence for one interesting parameter, Ax 2 = 2.706) is higher 
than that of the ISM in the LMC (Russell & Dopita 1992). There is no improvement (at greater 
than 90% confidence) in the fit when we allow the abundance of the Fe-like elements to vary. 

Although the model spectrum can still be rejected statistically at ~ 99% confidence, we con- 
sider the fit to be reasonably good. The model is certainly still an over-simplification. For example, 
we have assumed that the normalization of the soft bubble component is the same for both the east- 
ern and western parts of 30 Dor C. This assumption is problematic, however, since the power- law 
components of the two parts are very different. Because of the heavy absorption, the data below 
1 keV does not allow us to independently constrain the bubble component in the western part of 
30 Dor C. Dividing the data into the two parts (east and west) is an oversimplification; in general, 
the transition from east to west must be gradual and so the parameters inferred from our model 
should be considered as averaged values within the divided regions. We have also not included 
any systematic error in the data due to uncertainties in the calibration, background subtraction, 
and modeling, which are difficult to quantify at present. We therefore decided not to add more 
complications (such as varying abundances of individual elements) into the model. 



4. PHYSICAL PROPERTIES 

The morphology of the diffuse X-ray emission in 30 Dor C is consistent with that expected 
from the evaporation of cool gas into the hot interior of a bubble created in the ISM by the winds 
of massive stars. The close confinement of the limb-brightened X-ray emission by the outer Ha- 
emitting shell indicates that the bubble is roughly spherical. Otherwise, the bubble has to be 
elongated almost exactly along the line of sight. We have therefore interpreted the X-ray emission 
in terms of the stellar wind-blown bubble model of Weaver et al. (1977). 

Assuming that the above wind-blown bubble model provides a reasonable characterization 
of 30 Dor C, the physical properties of the enclosed hot gas can be inferred from the spectral 
parameters obtained above. We have adopted a notation similar to Shull &: Saken (1995), where 
te = t/ (10 6 yr) is the age of the superbubble in Myrs, no is the ambient density of the ISM in cm~ 3 , 
£38 = -^mech/ (10 38 erg s" 1 ) is the equivalent mechanical power injected into the superbubble in units 
of 10 38 erg s _1 , and kq is the ratio of the coefficient of thermal conduction to the classical Spitzer 
value (k sp = 6x 10~ 7 T 5 / 2 erg s _1 K _1 cm -1 ). The respective central temperature and density of the 
hot gas are T c = (5.3 x 10 6 K){L| 8 ng/(t|4°)} 1 /35 and Uc = ( 16 x 10 -2 cm -3^ L 6 8Tl i9 /c io/ i 22j.i/35 j 

assuming a superbubble radius R = (66pc){-L3stg/no} 1//5 (Shull & Saken 1995; see also Mac Low 
& McCray 1988). We can rearrange the above equations for T c , n c , and R to give 

(n c /1.6 x 10- 2 cm -3 ) (T c /5.3 x 10 6 K) +2 
n ° " (i?/66pc) 2 * 6 ' (3) 
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(n c /1.6 x 10~ 2 



cm" 3 ) (i?/66pc) 3 (r c /5.3 x 10 6 



K) 



(4) 



and 



Kq = 



(tc c /1.6 x l(T 2 cm- 3 ) (fi/66pc) 2 
(T c /5.3 x 10 6 K)5/2 t6 



(5) 



The parameters on the right side of equations (3)— (5) can be inferred from either the spectral 
fits with the bubble model (n c and T c ) or other independent measurements (R and ig). The central 
temperature and normalization of the best-fit bubble model (model vB+P) are T c = 7.4 x 10 6 K and 
K = 6.1 x 10 10 cm" 5 , respectively. This value of K gives a central density of n c = 1.4 x 1CP 2 cm -3 
and a mean emission measure of ~ 0.08 cm" 6 pc, assuming the bubble to be a prolate spheroid of 
dimension 50 pc x 45 pc. Adopting a mean radius of 47 pc for the bubble, which is the radius of the 
sphere that has the same volume as the prolate spheroidal bubble, gives an ambient density of n$ = 
(2.4 cm" 3 ) t|, a mechanical luminosity of L mec h = (4.3 x 10 erg s _1 )/tg, and a coefficient of thermal 
conduction of kq ~ 0.19 /tg of the Spitzer value. The mass of gas that has been swept up is then 
Mg = finoV ~ (3.3 x 10 4 Mq) tg, where fx = (14/11) tor is the mass per particle, mi is the mass of a 
hydrogen atom, and M Q is the mass of the Sun. Following Shull & Saken (1995), we also estimate 
the mass of X-ray emitting gas at temperatures ^ 10 5 K as Mx — (375/156)/in c V ~ 420 Mq. 
Finally, the thermal and kinetic energy of the superbubble are £th — (5/11) L mec h^ — 6.2 x 10 50 erg 
and Eke — (15/77) L mcc h£ — 2.7 x 10 50 erg respectively, and the thermal pressure of the hot gas 
is Pth — (2/3) (Eth/V) — 3.3 x 10 -11 dyne cm -2 . From the total energy, we infer that one or 
two core-collapse supernovae could have occurred within the superbubble, since each supernova 
supplies roughly 10 51 erg (e.g., Woosley & Weaver 1986). W-R stars are the most probable core- 
collapse supernovae progenitors in LH 90, and we may expect these stars to have a mass distribution 
similar to that of the W-R stars in the Galaxy, i.e., masses in the range from ~ 5 to ~ 30 Mq (e.g., 
Cherepashchuk 1991). Stars with masses in the range from ~ 12 to ~ 30 Mq are expected to 
undergo nuclear fusion of the a-elements and yield between ~ 0.15 and ~ 4 Mq of oxygen (Woosley 
& Weaver 1995; Thielemann et al. 1996). From the best-fit a-process elemental abundance of 1.6 
solar, we infer an oxygen/hydrogen mass ratio of 0.022 for the X-ray emitting gas in the bubble. 
The ISM in the LMC has an a-process elemental abundance of ~ 0.5 solar or an oxygen/hydrogen 
mass ratio of 6.7 x 10~ 3 . Thus, the bubble must have been enriched by ~ 5 Mq of oxygen, assuming 
that ~ 75 percent of the mass of the X-ray emitting gas is hydrogen. At least 2-3 high (^ 20 Mq) 
mass, core-collapse supernovae are needed to explain the oxygen abundance in the bubble, and the 
total energy supplied exceeds the value inferred from the spectral fits. 



5. CHECKING THE PHYSICAL PARAMETERS WITH MULTIWAVELENGTH 

DATA 



To make our results a bit more quantitative, we need to know the age of the bubble, which 
is still uncertain. Following Chu et al. (1995a), we use the observed expansion velocity of V cxp = 
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(0.59kms _1 ) R pc /t6 ~ 45 km s _1 (Dunne et al. 2001), where R pc ~ 47 in units of pc, to derive 
a bubble age of t ~ 6.2 x 10 5 years. This value of t is probably a lower limit to the true age of 
the bubble since the expansion is only detected in parts of the bubble (Y.-H. Chu 2004, private 
communication). The energy input from the OB association also tends to increase with time (Shull 
& Saken 1995). The ages of the most massive stars in the OB association give an age of ~ 3-4 Myr 
(Testor et al. 1993), although there are subclusters that are older (up to ~ 7 Myr). This kind of 
progressive star formation is also seen in the central OB association of the nearby 30 Dor nebula, 
typically with increasing strength; i.e., the most massive stars are formed later (e.g., Selman et al. 
1999). Therefore, we adopt a mean age of the OB association as ~ 4 Myrs. The inferred ambient 
density of no = 38cm~ 3 is reasonable for the expected high density ISM near 30 Dor C. 

The mean emission measure of the Ha-emitting shell is ~ 2000 cm~ 6 pc (Dunne et al. 2001). 
Adopting a minimum path length of 20 pc through the Ha emission region, which is twice as large 
as the width of the optical line-emitting filaments in Figure 2 (see also Fig. 11 of Dunne et al. 
2001), gives an upper limit to the mean electron density of n c ^ 10^0^ cm ~ 3 f° r the ionized shell, 
where £hq is the filling factor of the Ha-emitting gas. The thermal pressure of the ionized shell is 
approximately 2n e kT c ^ 2.8 x 10 ^ £h<!/ 2 dyne cm~ 2 , where T c ~ 10 4 K and k is the Boltzmann 
constant. This value is within a factor of three of the thermal pressure derived from the X-ray data 
for filling factors in the range 0.1-1. 

The mechanical luminosity required to power the bubble is -L mcc h = 1-1 x 10 ergs" 1 . There 
are 26 spectroscopically indentified O-type stars and 7 confirmed W-R stars in LH 90 (Testor et 
al. 1993). While the mass loss rates of W-R stars can be an order of magnitude higher than 
those of O-type stars (e.g., Nugis & Lamers 2000), the lifetime of these stars are fractions of a 
million years, i.e., much shorter than the probable age of the bubble, and so their energy input 
to the bubble is negligible compared with that from the O-type stars. Adopting a mean mass 
loss rate of M ~ lO^Moyr" 1 (de Jager, Nieuwenhuijzen, & van der Hucht 1988) and terminal 
velocities in the range from ~ 1000 to 3000 km s" 1 (Prinja, Barlow, & Howarth 1990), which 
are typical values for the O-type stars in LH 90, we calculate an integrated wind luminosity of 
L w = £(1/2)MV^ ~ 1-7 x 10 37 ergs _1 . Therefore, the O-type stars in the stellar clusters of 
LH 90 could easily supply the mechanical luminosity required to power the bubble, without the 
contributions from a few supernovae. If the bubble is much older than 4 Myr, then the mechanical 
energy input to the bubble from the O-type stars could greatly exceed the value of -L mec h inferred 
from the spectral fits. The best-fit bubble parameters give an estimate of the total energy input of 
-^mech* — 1.4 x 10 51 erg. This energy input is inferred from the density and temperature structure 
of the bubble, and not from its mechanical luminosity and age. 

There is clearly discrepancy between the expected mechanical energy input from the OB asso- 
ciation and that inferred from the bubble model. This apparent oversupply of mechanical energy is 
a common problem in the existing studies of several other structures around OB associations (e.g., 
Oey 1996; Naze et al. 2001; Cooper et al. 2003). In all these studies, including the present one, 
a potentially very important energy mechanism — dust grain cooling — is not considered. One 
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expects dust grains to be mixed with the X-ray-emitting gas due to the thermal evaporation and 
supernova ejecta. With a normal interstellar dust grain population, the cooling rate of dust grains 
heated by collisions with hot electrons can be substantially greater than the radiative cooling of 
the hot gas itself (e.g., Dwek k Arendt 1992 and references therein). Of course, dust grains subject 
to the sputtering. Their survival timescale depend on both the grain sizes and the ambient hot gas 
density. With the large uncertainties in these parameters, we find that it is conceivable for dust 
grains to be the major coolant of a superbubble. 

We find a coefficient of thermal conduction of kq = 0.05 ± 0.04 of the classical Spitzer value. 
Thermal conduction can be reduced by the presence of magnetic fields (e.g., Chandran k Cowley 
1998): the thermal conductivity can approach ~ 20% of the Spitzer value due to the chaotic 
transverse wandering of the magnetic field lines, although this value could be uncertain by a factor 
of two (Narayan k Medvedev 2001). The theoretical value of kq is consistent with our measurement, 
given the uncertainties in the bubble age. Thermal conduction is unlikely to be saturated (e.g., 
Cowie k McKee 1977), since the mean free path of a charge particle is much smaller than the radius 
of the bubble. 

As expected, the best-fit equivalent hydrogen column densities toward the eastern and western 
parts of 30 Dor C are significantly higher than the equivalent Galactic H I column density toward 
30 Dor C of A H (Gal) = 6.4 x 10 20 cm" 2 (Dickey k Lockman 1990). We therefore use the optical 
extinction of stars in the OB association to independently estimate the total (atomic plus molecular) 
column density toward 30 Dor C. The mean interstellar extinction for the association as a whole 
is Eb-v ~ 0.4 mag and reaches a maximum of E-q-y ~ 0.8 mag in the western part (Testor et 
al. 1993). Adopting a value of E-q-v = 0.08 mag for the Galactic foreground extinction toward 
30 Dor C (Xu et al. 1992) and assuming gas-to-dust ratios of Ah/£b-v = 4.8 x 10 21 cm" 2 mag" 1 
(Bohlin, Savage, k Drake 1978) and N H /E B -v = 2.4 x 10 22 cm" 2 mag" 1 (Fitzpatrick 1985) for 
the Galaxy and LMC respectively, we estimate a mean column density of Ah ~ 8.1 x 10 21 cm" 2 
toward the OB association and a maximum column density of Ah ~ 1.8 x 10 22 cm" 2 toward the 
western part. While these values appear to be larger than our measurements, we had assumed 
solar abundances for the atomic absorption cross-sections in the spectral fits. Almost all of the 
absorption cross-section at ~ 1 keV is provided by metals and so, for a metal abundance of 0.5 
solar, the measured X-ray absorption column densities towards the eastern and western parts of 
the bubble are Ah ~ 10 22 cm" 2 and ~ 2 x 10 22 cm -2 , respectively. These X-ray absorption column 
densities are consistent with the absorption column densities estimated from the optical extinction 
of stars in the OB association. The difference in the measured X-ray absorption column densities of 
Ah ~ 10 22 cm" 2 is consistent with the H2 column density inferred from the CO emission projected 
onto the western part of 30 Dor C (see § 3.1). 
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6. NATURE OF THE NONTHERMAL X-RAY EMISSION 

6.1. Synchrotron Emission 

For energy-conserving shocks, the Ha expansion velocity of 45 km s _1 (Dunne et al. 2001) 
corresponds to a shock velocity of ~ 60 km s~ x (e.g., Chu et al. 1995b), and gas passing through 
the shock would only be heated to temperatures of ~ 10 5 K (e.g., Hollenbach & McKee 1979), 
which is well below the temperature of the X-ray emitting gas. The observed [S ll]/Ha flux ratio 
of 0.3 (Mathewson et al. 1985) is also below that expected from shocked gas (e.g., Long et al. 
1990) and, contrary to the observations, the bulk of the nonthermal X-ray emission is expected to 
originate in the eastern part of 30 Dor C since the non-thermal radio emission is much stronger 
on the eastern side of the bubble than on the western side. Thus, we conclude that any recent 
supernova explosion must have occured deep inside the bubble, producing little or no observable 
nonthermal X-ray synchrotron radiation. Bamba et al. (2003, 2004) explained the the apparently 
nonthermal component in the X-ray spectrum of 30 Dor C as synchrotron radiation from relativistic 
electrons. If this is true, the required electron energies should be a few hundred TeV for a typical 
interstellar magnetic field of order B = 3 x 10~ 6 G. Although low (up to a few tens of km s _1 ) 
velocity shoulders are visible on the Ha emission lines from various parts of the bubble (Y.-H. Chu 
2004, private communication), there is a complete lack of emission with velocities greater than 100 
km s _1 from the recessional velocity of the LMC and so the shock velocities are certainly too low 
to produce such energetic particles. 

The most obvious site for an effective particle acceleration should be the reverse shock, which 
separates the free-streaming stellar wind from the region of thermally evaporated material (R\ in 
Fig. 1 of Weaver et al. 1977). For a mean mass loss rate of M ~ 10~ 6 M Q yr _1 , a terminal velocity 
of Voo ~ 2000 km s -1 , an ambient density of no = 38 cm~ 3 , and a bubble age of 4 Myrs, the inner 
shock radius is roughly 7 pc (Weaver et al. 1977). The time scales for the relativistic electrons 
to diffuse from the reverse shock to the edge of the bubble would have to be shorter than their 
synchrotron half-lives, which are of order (5 x 10 3 years) (B/3 x 10~ 6 G)~ 3 / 2 (e.g., Lang 1974, pp. 
31-32). The relativistic electrons would therefore have to travel in the bubble at velocities that 
are an order of magnitude faster than the sound velocity of ~ 10 2 -10 3 km s _1 and two orders 
of magnitude faster than the Alfven velocity of S/(9.67run c ) 1 / 2 ~ 30 km s _1 , which is untenable 
(e.g., Jaffe 1977). It is possible that the electrons could be accelerated by turbulent motions in 
the bubble itself. The turbulent velocities would be of order or greater than the Alfven velocity, 
since otherwise the magnetic field would damp the turbulent motions (e.g., Jaffe 1977), and so the 
turbulent kinetic energy involved would be a few 10 48 erg. The electrons must be accelerated on 
time scales that are shorter than their radiative lifetimes and so the power dissipated in this process 
is of order 10 37 erg s _1 . While this value is comparable to the mechanical wind luminosity of the 
early- type stars, the "second-order" acceleration by turbulence is probably too inefficient. 
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6.2. Inverse Compton Radiation 



The inverse Compton (IC) scattering of Cosmic Microwave Background (CMB) and other 
ambient photons (e.g., infrared photons from reprocessed starlight) by relativistic electrons will 
also generate nonthermal X-ray emission. The upscattering of such photons to 5 keV requires 
electrons with energies in the range from ~ 0.1 to ~ 1 GeV. These electrons radiate synchrotron 
emission at frequencies of MHz to tens of MHz in mag netic fields of order 3 x 10" 6 G. For CMB 
photons, the ratio of inverse Compton-scattered X-ray flux to radio synchrotron flux is 



(e.g., Blumenthal & Gould 1970), where S x and S T are the flux densities at frequencies v x and v T , 
respectively. The functions b(s) and a(s) are given in Blumenthal & Gould (1970) for a power-law 
distribution of electron energies with index s. The radio synchrotron spectrum should have the same 
spectral index as the IC emission and so we have assumed a spectral index of a = (s — l)/2 = 1.5 
for the radio synchrotron spectrum. The ratio of b(s) to a(s) is ~ 420 for this value of a. We 
have also assumed that the magnetic field strength is 3 x 10~ 6 G and that the electron distribution 
is an unbroken power-law extending from Lorentz factors of few times 10 3 , which will upscatter 
CMB photons to X-ray wavelengths, to Lorentz factors of order 10 4 , which will produce 843 MHz 
radio emission via synchrotron radiation. This last assumption may not be correct, however, since 
the nonthermal X-ray emission is stongest in the western part of the bubble where there is no 
evidence for nonthermal radio emission (Mathewson et al. 1985). Nonetheless, the total (thermal 
plus synchrotron) emission at 843 MHz is 1.2 Jy (Mathewson et al. 1985) and so the expected 1 keV 
to 843 MHz flux density ratio of ~ 1.4 x 10~ 6 (B/3 x 10~ 6 G)~ 5/2 gives an upper limit to the IC 
emission at 1 keV of ~ 2.5 x 10~ 3 (B/3 x 10~ 6 G)~ 5 / 2 photons cm 2 s 1 keV 1 . Although this 
value is consistent with the flux densities of the power-law emissions from the eastern and western 
parts of 30 Dor C, the IC emission is probably much fainter than this, since the flat (a ~ 0) radio 
spectrum is an indication that the radio emission is mostly thermal (Mathewson et al. 1985). 

The upscattering of infrared photons to X-ray energies requires electrons with Lorentz factors 
of a few hundred. For a power-law distribution of electron energies with index s ~ 4, we would 
expect the number of electrons capable of upscattering infrared photons to X-ray energies to be 
~ 10 2 times the number of electrons capable of upscattering CMB photons to X-ray energies. The 
42.5-122.5 fim luminosity of 30 Dor C can be estimated from 



(Cataloged Galaxies and Quasars Observed in the IRAS Survey, Appendix B), where feo and /ioo 
are the IRAS flux densities in Jy at 60um and 100/iin, respectively, and D is the distance to the 
source in cm. The extended infrared emission measured by IRAS is roughly 15 and 30 Jy at 60 and 
100/itm, respectively (Laspias k, Meaburn 1991), and so we estimate a 42.5-122.5 fim luminosity of 
the order of a few 10 38 erg s _1 . Thus, the mean infrared energy density of the bubble is of the order 




(6) 



L fir = AirD 2 1.26 x 10- 11 (2.58/ 60 + /ioo) ergs' 



-i 



(7) 
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10~ 13 erg cm -3 and could be higher if the infrared emission is concentrated in the outer shell. If 
we adopt a mean photon energy of e = 0.018 eV for the infrared emission, then the total infrared 
photon number density is ~ 3.5 cm~ 3 . The CMB photons have, in contrast, a mean photon energy 
of e = 6.25 x 10~ 4 eV and a total photon number density of 400 cm' 3 (e.g., Gaiser, Protheroe, & 
Stanev 1998). From equation (6) of Gaisser et al. (1998) and assuming that the electron energies 
are distributed as a power-law of index s ~ 4, we find that the IC emission at a fixed photon 
energy scales roughly as e 3 / 2 n p h, where n p h is the photon number density of the radiation under 
consideration. So we would expect the X-ray emission from IC-scattered infrared photons to be 
~ 30 percent higher than that from the upscattering of the CMB photons. However, the optical 
radiation from stars should have an energy density greater than that of the infrared emission and 
could be IC-scattered to X-ray energies by electrons with Lorentz factors of ~ 50. If we assume 
that the mean energy of the optical photons is e ~ 3.1 eV, then the total photon number density 
of the optical radiation should be ^ 0.02 cm" 3 and the X-ray emission from IC-scattered optical 
photons should be an order of magnitude higher than that from the upscattered infrared and CMB 
photons. Unless the particle energy distribution is far from the power law (e.g., a Maxwellian 
distribution) and/or the particle population is greatly enhanced at the shell, the nonthermal X-ray 
emission should then be concentrated in the central region of 30 Dor C, where the optical emission 
is greatest. This is inconsistent with the observed limb-brightened hard X-ray emission (Fig. 3). A 
strong shock, as expected inside the bubble, will also produce electron spectra with indexes of s ~ 2 
(e.g., Blandford & Eichler 1987), which is totally at odds with the spectral index of the nonthermal 
X-ray emission. On the other hand, if the relativistic particles represent the thermalized wind 
material from young pulsars during the lifetime of the bubble, then the particle energies could be 
distributed as a Maxwellian of peak energy 7 ~ 10 2 -10 3 (Arons & Tavani 1994; Atoyan 1999). 
Therefore, we cannot rule out the possibility that the bulk of the nonthermal X-ray emission from 
the bubble is due to IC radiation. 



6.3. Nonthermal Bremsstrahlung Radiation 

Nonthermal particles generate bremsstrahlung in regions where the energy losses due to ion- 
ization and Coulomb collisions are significant. The energy flux of the bremsstrahlung radiation is 
^ 10 -5 times the total energy input to the particles (e.g., Dogiel et al. 2002), i.e., a large amount 
of energy is required to maintain the particles against Coulomb/ionization losses to the emitting 
region. The best-fit bubble model (model vB+P) has power-law components with unabsorbed 0.1- 
10 keV luminosities of 1.2 x 10 36 and 2.8 x 10 36 erg s _1 in the eastern and western parts of the 
bubble, respectively. The required energy input is therefore <; 10 41 -10 42 erg s~ x . Although the 
necessary energy input can be less than this for a quasi-thermal distribution of electrons formed in 
regions of in situ acceleration (Dogiel et al. 2002), it is still orders of magnitude greater than the 
integrated wind luminosity of ~ 1-7 x 10 37 erg s _1 from the O-type stars in LH 90. A recent (within 
the last a few ~ 10 3 years; Chu & MacLow 1990) supernova explosion could provide the energy 
input to the bubble, but this would mean that we are seeing the event by chance. The 26 O-type 
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stars in LH 90 have masses in the range from 20 to 40 M (e.g., Walborn et al. 1995, 1999). If we 
assume a Spitzer initial mass function for the stars in LH 90, then we would expect ~ 12 O-type 
stars with masses in the range from 40 to 60 Mq. Although some of these stars would undoubtedly 
have evolved into W-R stars, the remainder would have exploded as supernovae (stars with masses 
greater than ~ 60 M Q do not explode as supernovae, but instead collapse into a black-hole). Thus, 
there must have been at least ~ 5-6 supernovae in 30 Dor C, and so the probability of seeing a 
supernova within the last 10 4 years is approximately 1 percent, assuming that the age of the bubble 
is 4 Myr. 

7. CONCLUSIONS 

In summary, the nature of the apparent nonthermal hard X-ray component is still uncertain. 
The expansion of the bubble is much too slow in order to produce the energetic particles necessary 
for nonthermal X-ray synchrotron emission. Although such energetic particles can be produced 
at the inner reverse shock, they cannot cross the bubble within their synchrotron half-lives. Non- 
thermal bremsstrahlung is an inefficient process, requiring at least three orders of magnitude more 
energy than can be supplied by the O-type stars in LH 90. A recent supernova could supply the 
necessary energy, but this would mean that we are seeing the event by chance (a probability of < 1 
percent). The IC scattering of CMB and infrared photons cannot account for all of the nonthermal 
X-ray emission unless the radio flux is purely synchrotron radiation, which seems unlikely. While 
the X-ray emission from IC-scattered optical photons should be much higher than that from IC- 
scattered CMB and infrared photons, it should be concentrated in the central region of 30 Dor C, 
where the optical photon density is greatest, in sharp contrast to the data. The nonthermal X-ray 
emission could instead be due to the IC scattering of CMB and infrared photons by relativistic 
particles originating in young pulsar wind material with energies peaking between 7 ~ 10 2 -10 3 . 
High spatial resolution X-ray imaging, which can be afforded by an on-axis Chandra observation, 
will be particularly useful to further the understanding of the component. Such an observation 
will allow for an accurate determination of the size and geometry of the X-ray-emitting region and 
for a detailed comparison with images in other wavelength bands. A substantial fraction of the 
expected energy input from stellar winds and supernovae may be radiated in the infrared by dust 
grains entrained in the hot gas. 

We thank R. Chris Smith for providing his optical image of 30 Dor C in computer-readable 
format, Glenn Allen for discussions on the nonthermal emission, and the referee for insightful 
comments on an earlier draft of the paper. This research is supported by the NASA LTSA grant 
NAG5-8935. 
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Table 1. Spectral Fits to the Integrated Emission from 30 Dor C in the 0.3-10 and 1-10 keV 

bands. 



Model* 







Parameter 


B+P 


M+P 


vB+P 


A^cast 15 




(10 21 cm" 2 ) 


(5 4+1-1 


c Q4+0.55 


4.0818:18 


^H.wcst 




(10 21 cm' 2 ) 


11 41+0.96 
il - 41 -0.71 


11.90lg;» 


io-i8l8:S 


kT d 




(keV) 


U.Od»_Q Qgg 


o-i86±g : gsi 


0.6418;?° 


Z e 




(Ze) 


0.50 f 


0.50 f 


0.50 f 


Z s 




(Ze) 


0.30 f 


0.50 f 


i.63i8:g 


K h 
^thermal 




/ n e n H dV) 


(3.71?- 2 ) x 10- 3 


(1.5618:%) x io- 2 


(3.1+ 3 ; 2 ) x 10- 4 


r 




2 79+ - 11 
L - 1 y -0.10 


o oc+0.10 

z - 00 -o.io 


z - 01 -0.11 


^pl,east^ 


(io- 


" 4 photons cm -2 s -1 keV -1 ) 


°- yo -0.54 


6 41+ - 54 


4.3418:" 


K k 
Is - pl,west 


(io- 


~ 4 photons cm -2 s -1 keV -1 ) 




13.5li3 


9.611:1 


X 2 (dof) 






459.8(347) 


463.7(347) 


430.6(346) 


Prob. 1 






4.5 x 1(T 5 


2.7 x 10~ 5 


1.3 x 10~ 3 
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a Model: B = wind-blown bubble (model valid for gas temperatures in the range from 8 eV to 
80 keV); vB = wind-blown bubble with variable abundances (model valid for gas temperatures in 
the range from 8 eV to 80 keV); M = MEKAL (model valid for gas temperatures in the range from 
8 eV to 80 keV); P = power-law. 

b Measured column density towards the eastern part of the bubble. 

c Measured column density towards the western part of the bubble. 

d The gas temperature at the center of the bubble. 

C C, N, Na, Al, Fe, and Ni abundance relative to solar of the bubble. 

f Fixed at the known value of 30% solar for the ISM in the LMC. 

g O, Ne, Mg, Si, S, Ar, and Ca abundance relative to solar of the bubble. 

h The normalization of the bubble model. The units are cgs. 

'Photon spectral index. 

J Normalization for the power-law model in the eastern part of the bubble. 
k Normalization for the power-law model in the western part of the bubble. 

Probability that the model describes the data and that \ 2 exceeds the observed value by chance. 
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RIGHT ASCENSION (J2000) 

Fig. 1— The XMM-Newton image of the 30 Dor C field in the 0.3-8 keV band with the bright 
sources of X-ray emission labelled. The image has a pixel size of 3" x 3" (0.72 pc x 0.72 pc). The 
solid black lines mark source and background (Bkgd 1 for the first observation; Bkgd 2 for the 
second observation) regions from which spectra were extracted. The shading is proportional to the 
logarithm of the intensity and ranges from 10 -3 (white) to 830 (black) counts pixel -1 . 
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RIGHT ASCENSION (J2000) 



Fig. 2. — An adaptively smoothed XMM-Newton image (contours) of 30 Dor C in the 0.3-8 keV 
band superposed on an Ha MCELS image (grayscale) at 6563A (Smith et al. 1998). In the adaptive 
smoothing process, exposure corrected images are convolved with a two-dimensional Gaussian, with 
its width adjusted to achieve a count-to-noise ratio greater than 6. The shading is proportional to 
the logarithm of the intensity and ranges from 57 (white) counts pixel -1 to 4200 (black) counts 
pixel -1 in the optical image. Countours are drawn at (1.4,1.9,3.0,4.8,8.9, and 33) x 10 -2 counts 
s -1 arcmin -2 in the X-ray image. 
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RIGHT ASCENSION (J3000) RIGHT ASCENSION (J3000) RIGHT ASCENSION (J3000) 



Fig. 3. — Adaptively smoothed XMM-Newton images of the central region of 30 Dor C field in 
the (a) 0.3-1, (6) 1-2, and (c) 2-8 keV bands. The X-ray images have been adaptively smoothed 
using the same procedure as in Fig. 2. The shading is proportional to the logarithm of intensity 
and ranges from 3.2 x 10~ 3 (white) to 0.48 (black) counts s -1 arcmin -2 in the 0.3-1 keV, from 
2.7 x 10~ 3 (white) to 0.41 (black) counts s -1 arcmin -2 in the 1-2 keV, and from 3.1 x 10~ 3 (white) 
to 0.14 (black) counts s _1 arcmin~ 2 in the 2-8 keV band images. Contours are drawn at (4.8, 6.4, 
10, 16, 30, and 110) xl0~ 3 counts s" 1 arcmkr 2 in the 0.3-1 keV, at (4.1, 5.4, 8.6, 14, 25, and 93) 
xlO -3 counts s _1 arcmin^ 2 in the 1-2 keV, and at (4.6, 6.1, 9.7, 15, 28, and 100) xl0~ 3 counts 
s _1 arcmin -2 in the 2-8 keV band images. 
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Fig. 4. — The co-added MOS spectra of the 30 Dor C superbubble. The upper panel shows the data 
together with the folded model (solid lines) comprising the best-fit wind-blown bubble spectrum 
(dashed lines) and power-law continuum (dotted lines). The residuals, in units of a, from the best- 
fit model are shown in the lower panel. The data have been binned so that the signal-to-noise ratio 
in each bin exceeds 4. (a) Eastern part of superbubble. (b) Western part of superbubble. 



